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Abstract. Very massive primordial stars (140 Mq < M < 260 Mq) are supposed to end their 
lives as PISN. Such an event can be traced by a typical chemical signature in low metallicity 
stars, but at the present time, this signature is lacking in the extremely metal-poor stars we are 
able to observe. Does it mean that those very massive objects were not formed, contrarily to the 
primordial star formation scenarios ? Could it be possible that they avoided this tragical fate ? 

We explore the effects of rotation, anisotropical mass loss and magnetic field on the core size 
of very massive Population III models. We find that magnetic fields provide the strong coupling 
that is lacking in standard evolution metal-free models and our 150 Mq Population III model 
avoids indeed the pair-instability explosion. 
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1. Introduction 

According to lHeger et al. (2003), the fate of single stars depends on their He-core mass 



(M a ) at the end of the evolution. They have shown that at very low metallicity, the stars 
having 64 Mq < M a < 133 Mq will undergo pair-instability and be entirely disrupted 
by the subsequent supernova. This mass range in M a has been related to the initial 
mass the star must have on the main sequence (MS) through standard evolution models: 
140 Mq < Mi n i < 260 Mq. Since we will present here a non-standard evolution, we will 
rather keep in mind the M a range. 

The typical mass of Population III (Pop III) stars is explored by early struc ture forma- 



tion studies and chemistr y considerations about cooling. Different studies (see I Abel et al 



20021 : Bromm et al.ll2002l among others) give the same conclusion: Pop III stars are sup- 



posed to be massive or very massive, eve n when a bimodal mass distribu tion allows the 
formation of lower mass components (see iNakamura and Umemura 200 ll ). Therefore we 



expect that many among them should die as pair-instability supernovae (PISN) . 

1.1. A typical chemical signature which remains unobserved 
These PISN events are su pposed to leave a typical chemical signature. According to 



Heger and Wooslevl (|2002l ). the complete disruption of the star leads to a very strong 
odd-even effect: the absence of stable post-He burning stages deprives the star of the 
neutron excess needed to produce significant amounts of odd-Z nuclei. Also the lack of 
r- and s-process stops the nucleosynthesis around zinc. Even if one mixes these yields 
with the yields of zero-metallicity 12 — 40 Mq models (which end up as standard Type 
II SNe), the PISN signature re mains and should be ob servable. 

However, using those yields. rTumlinson et al. ( 2004f ) have shown that it provides only 



a very poor fit to the abundances pattern observed in the metal-poor stars known today. 
The odd-even effect is not observed, and the models significantly over-produce Cr and 
under-produce V, Co and Zn. 
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The most metal-deficient stars are supposed to be formed in a medium enriched by 
only one or a few SNe. The absence of the chemical signature of the PISN is a strong 
argument against their existence. But how could that be? 

1.2. Simple solutions 

The simplest solution to explain this absence is to suppose that the mass domain in 
question was not formed in the primordial clouds. Maybe the primordial IMF was not as 
top-heavy as we actually think, and the most massive stars formed then could very well 
be too small for such a fate. H owever, recent works on primordial stars form ation seem 
to rule out this possibility (see Yoshida et al. 20061 : lO'Shea and Normanll2007 ). 



Another possibility is that the signature was very quickly erased by the next gener- 
ations of stars. Maybe the metal-poor stars we observe are enriched by more SNe than 
we actually think, and the later contributions are masking the primordial ones. Only 
the observation of more and more metal-deficient stars will provide an answer to that 
possibility. 

One can also wonder whether there would be a way for those stars to avoid their fate. 
In this context, the simplest solution is to suppose that some mechanism could lead to 
such a high mass loss that the conditions for pair-instability would no more be met in 
the central regions. The aim of the present work is to explore this possibility. 



2. Rotation at very low metallicity 

Radiatively-driven winds are supposed to scale with metallicity as ^(°- 5 -°- 86 ). Thus, 
very low or even no mass lo ss is expected when Z approaches to 0. 

However, we have shown ( Mevnet et al. 2006h that rotation can change the mass loss 
history of low metallicity stars in a dramatic way. Two processes are involved: 

(a) During MS, because of the low radiative winds, the star loses very little mass and 
thus very little angular momentum. As the evolution proceeds, the stellar core contracts 
and spins up. If a coupling exists between the core and the envelope {i.e. meridional 
currents or magnetic fields), the surface may be accelerated up to the critical velocity 
and the star may lose mass by a mechanical wind due to centrifugal acceleration. The 
matter is launched into a decretion disk, which may be dissipated later by the radiation 
field of the star. 

(b) Rotation induces an internal mixing that enriches the surface in heavy elements. 
The effective surface metallicity is enhanced by a factor that has been shown to be very 
large (up to 10 6 for a 60 Mq at Z- m [ — 10~ 8 ). Rotation also favours a redward evolution 
after the MS, allowing the star to spend more time in the cooler part of the HR diagram. 





Z = 0.020 

Z-0.002 

z=io- 6 

Z = 



Figure 1. Left: Variations in the radius as a function of initial mass for various metallicities. 
Centre: Amplitude of the radial component of the meridional circulation inside 20 Mq models 
at the same evolutionary stage (X c = 0.4) and having the same rotation rate (v/v CI it = 0.5) but 
various me tallicities. Right : Evol ution of the equatorial velocity during the MS. The models are 
taken from lEkstrom et al.1 (|200Sfl . 
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The opacity of the envelope is increased, and the radiative winds may thus be drastically 
enhanced. 

Both effects add up and lead to strong mass loss at very low metallicity. 

2.1. Rotation at Z = 

But what happens when the metallicity drops down to Z = strictly? The absence of 
carbon prevents the star to start burning hydrogen through CNO-cycle, but pp— chains 
are not energetic enough to sustain the star, so it contracts longer during its formation. 
At metallicities Z = 0.020, 0.002, 10" 5 , and respectively, the radii on the ZAMS will 
be in a ratio 3.5 : 2.8 : 2 : 1 for a 20 Mq model (see Fig.[TJ left), that is the Z = star will 
be twice as compact as the Z = 10 -5 one. This has a direct influence on the amplitude 
of the outer cell of the meridional circulation, which depends on the compactness of the 



star as U r 



1 



12- 



For the same metallicities, the amplitude of the outer cell will 
4 : 1 (Fig. [TJ centre). This means that at Z = 0, the core- 



2-r Gp 

be in a ratio 100 : 17 
envelope coupling is almost null. Figure [T] (right) shows the evolution of the equatorial 
velocity during the MS for our 20 M© models. At standard metallicity, the meridional 
circulation is strong, but the mass loss is also strong, so the model spins down. At lower 
but non-zero metallicities, the mass loss drops, and though the core-envelope coupling is 
weak, the models spin up. At Z = 0, the mass loss is null, but the model evolves almost 
in a regime of local conservation of angular momentum. Thus, the natural inflation of 
the radius during the MS leads to a spin down of the model. The model may reach the 
critical limit, but later in its MS evolution (for a given initial mass and velocity), and 
the mechanical mass loss remains very small. 

At central H exhaustion, the core of the Z — model is already hot enough to burn 
some helium, so the transition between H-burning and He-burning is smooth, without any 
structural readjustment. The model remains in the blue part of the Hertzsprung-Russell 
diagram until very late in the central He-burning phase (Fig. [5]). When it eventually 
evolves towards the red part, the outer convective zone is very thin, so the enrichment of 
the surface remains very low. Whereas at Z = 10~ 8 we get a strong post-MS mass loss, 
at Z = it is negligible. 

Rotation alone seems thus to fail in providing a way to lose sufficient mass at Z — 0. 

2.2. Two natural effects of rotation 

In the previous calculations, we have neglected two mechanisms that arise naturally with 
rotation and could change this picture. 

The first one is the wind anisotropy (see Maeder 1999): when a star rotates fast, 



. Z=10-' 
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Figure 2. Evolution of T c ff during core He-burn ing for two mo dels of 60 Mq. The Z = 10 

model is taken from iHirschi ( 20071). 



4 



Ekstrom & al. 



it becomes oblate. As shown by von Zeipel theorem, the poles become hotter than the 
equator, and the radiative flux is no more spherically symmetric: it gets much stronger in 
the polar direction than in the equatorial plane. Since the mass is lost preferentially near 
the poles, it removes much less angular momentum than in the spherical configuration, 

and thus the star may reach the critical limit much earlier. 

The second mechanism is the magnetic fields. According to Spruit! ( 2002 ) the differen- 
tial rotation may amplify an existing magnetic field through the Tayler-Spruit dynamo 
mechanism, and provide a strong core-envelope coupling. This coupling will be able to 
accelerate the surface very early in the evolution, and maintain the star at critical limit 
throughout its entire evolution. 



3. Evolution with anisotropical winds and magnetic fields 

3.1. Ingredients 

With the help of the two effects mentioned above, we have computed an exploratory 
model. Its mass has been chosen to be 150 M & and the initial ratio between the equatorial 
velocity and the critical one to be Umi/w C rit = 0.56. Let us mention that this ratio is higher 
than the one needed to account for the observed average velocities at solar metallicity, 
but it is still a "reasonable" value, not an extreme one. 

The computation has been accomplished using the Geneva code with up to date nuclear 
reaction rates obtained with NETGEN (http://www-astro.ulb.ac.be/Netgen/). The 
opacity tables come from OPAL (http : // www-phys . llnl . gov/ Research/OPAL/opal .html) 
with the extension at low temperature bv lFerguson et al" (2005). The initial composition 
is X = 0.753, Y = 0.247 and of course Z = Q. 

The radiative mass loss prescri ption is an i mpor tant ingredient of the modelization of 
massive stars. Here we have used iKudritzkil (|2002l) . Since this pres cription is not aime d 
at the case Z = strictly, we have used the same ad aptations as inlMarigo et al. ( 20031 ). 
The Wolf-Rayet (WR) mass loss rate is taken fr om Nugis and Lamersl (|2000l )"~with the 
metallicity scaling from lEldridge and Vinkl ( 2006f ). For the calculation, we have taken the 
effective surface metallicity Z e g = (1 — X — Y) sur f so that the enrichment of the surface 
is accounted for. We must stress that this Z e g is mainly composed by CNO elements but 
no iron. It is usually considered that WR winds are triggered by Fe lines, whereas the 
CN O lines determine only v^ , so we expect that no WR winds can take place at Z = 0. 
But Vink and de Koterl ( 2005 ) have shown that when the metallicity gets really low, the 
CNO lines take over the role of Fe lines in the line driving, and this is what we have 
assumed here. 

When the m odel reaches the crit ical limit, the mechanical mass loss has been treated 
as described in lMevnet et al. | (|2006h . 

The treatment of anisotropic winds has been implemente d as in iMaederl ( 2002f ) and 
the effect of magnetic fields as in iMaeder and Mevne 3 (|2005l) . 



3.2. Result 

In Fig. [31 we present the evolution in the HR diagram (left panel) and the evolution of 
mass with time (right panel). The grey line shows a non-rotating model computed with 
the same physics for comparison. 

During its whole evolution up to the end of core He-burning, the non-rotating model 
loses only 1.37 Mq. This illustrates the weakness of radiative winds at Z = 0. 

The evolution of the rotating model (black line) can be described by four distinct 
stages: 
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(a) (continuous part, lower left corner) The model starts its evolution on the MS with 
only radiative winds, losing only a little more than 0.002 Mq. During this stage, it 
accelerates quickly, mainly because of the strong coupling exerted by the magnetic helds. 

(b) (dashed part) When its central content of hydrogen is still 0.58 in mass fraction, 
it reaches the critical velocity and starts losing mass by mechanical wind. It remains at 
the critical limit through the whole MS, but the mechanical wind removes only the most 
superficial layers that have become unbound, and less than 10% of the initial mass is 
lost at that stage (11.44 Mq). The model becomes also extremely luminous, and reaches 
the Eddington limit when 10% of hydrogen remains in the core. Precisely, it is the so- 
called fir-limit that is reached here. Due to the fast rotation, the maximum Eddington 
factor allowed is reduced; at the same time, because of the high luminosity, the critical 
velocity is reduced in com parison with the one derived from the classical Jl-limit (see 
Maeder and Mevnetl[2000l for details). 

(c) ( dotted part) The combustion of helium begins as soon as the hydrogen is exhausted 
in the core, then the radiative H-burning shell undergoes a CNO flash, setting the model 
on its redward journey. The model remains at the Sir-limit and loses a huge amount of 
mass. The strong magnetic coupling keeps bringing angular momentum to the surface and 
even the heavy mass loss is not able to let the model evolve away from the critical limit. 
The mass lost during that stage amounts to 53.46 M Q . When the model starts a blue 
hook in the HR diagram, its surface conditions become those of a WR star (X sur f < 0.4 
and T e g > lO'OOO K). The luminosity drops and takes the model away from the r-limit, 
marking the end of that stage. 

(d) (continuous part) The rest of the core He-burning is spent in the WR conditions. 
The mass loss is strong but less than in the previous stage: another 26.34 Mq are lost. 

At the end of core He-burning, the final mass of the model is only Mg n = 58 Mq, 
already below the minimum M a needed for PISN (M a ^ 64 Mq). Note that the con- 
traction of the core after helium exhaustion brings the model back to critical velocity, so 
this value for Mfi n must be considered as an upper limit. 




Log T e „ [°K] age [Myr] 

Figure 3. Black line: rotating model; continuous part: beginning of MS (X c — 0.753 down to 
0.58; dashed part: rest of the MS; dotted part: beginning of core He-burning phase (Y c — 1.00 
down to 0.96); continuous part: rest of the He-burning. Grey line: non-rotating model for com- 
parison. Left panel: evolution in the Hertzspmng-Russell diagram; Right panel: evolution of 
the mass of the model. The mass indicated is the mass lost at each stage, not a summation. 
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4. Summary 

The model we presented here is exploratory. We cannot draw general conclusions from 
it. But our model shows that heavy mass loss is possible even at Z = 0, and the answer 
to our title's question is: yes, under certain conditions, very massive stars can indeed 
avoid PISN. 

Some aspects need yet to be clarified. Is the WR mass loss rate we have used really 
valid? Can the CNO lines alone really drive a WR wind? In a more general perspective, 
we still lack a good mass loss recipe for the strict Z — case. A word of caution must 
also be cast on the inclusion of magnetic fields. The validity of the Tayler-Spruit dynamo 
is still under debate, and more work need to be done before we may confidently rely 
on results that have been obtained with the actual treatment. Moreover, there is not 
a clear consensus whether magnetic fields were present in the early Universe or not. In 
the Taylor-Spruit dynamo, the mechanism amplifies a pre-existing field, so if none were 
present at the start, we cannot use it. 

Anyway, the physics used in the present model is today's "state of the art" and it 
is interesting to study what can be achieved with it. Our result is encouraging, because 
the computation has been accomplished with reasonable assumptions: the initial rotation 
rate used here was fast but not extreme, and the mechanisms called upon (anisotropy of 
the winds and magnetic fields) are not exotic ones, but two natural effects which arise 
when one treats properly the case of rotation. 

After this first step, more work is needed. First, we have to refine the mass loss treat- 
ment at fir-limit and check if the mass loss stays as strong as here. If it does so, we 
have to check if our result are valid in the whole PISN mass domain. Higher mass mod- 
els should experience higher mass loss, but it is necessary to check if it would still be 
sufficient to help them avoid pair-instabilities. 
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